Fictitious fluxes in doped antiferromagnets 
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In a tight binding model of charged spin-^ electrons on a square lattice, a fully polarized ferro- 
magnetic spin configuration generates an apparent (7(1) flux given by 27r times the skyrmion charge 
density of the ferromagnetic order parameter. We show here that for an antiferromagnet, there are 
two "fictitious" magnetic fields, one staggered and one unstaggered. The staggered topological flux 
per unit cell can be varied between — tt < $ < tt with a negligible change in the value of the effec- 
tive nearest neighbor coupling constant whereas the magnitude of the unstaggered flux is strongly 
coupled to the magnitude of the second neighbor effective coupling. 



GENERAL PROPERTIES OF SPIN GENERATED 
FLUXES 

It is well known that in the theory of the quantum HaU 
problem spin textures can generate an effective U{1) flux 
which acts as an effective magnetic field leading to an. as- 
sociation between topological and electrical charge. El In 
this paper we investigate to what extent this effect can 
be generalized to the antiferromagnet. 

For definiteness, we consider the t-J model on a square 
lattice. This model is described by the Hamiltonian 



{f^clc,,+h.c.) + J^\S,-S, 



1 



(1) 



where the summation runs over nearest neighbor pairs 
(ij) and — cl^ffa/iCifs and a denotes the vector of 
Pauli matrices. AU states containing doubly occupied 
sites have been excluded from the Hilbert space. [For a 
general reference, see Ref. ||.] 

To make explicit the connection between spin rotations 
and effective couplings, we introduce a local change of 
spin coordinates [see for instance Ref. [sj, choosing the 
local spin quantization axis at site i along fli. The ac- 
tion of this local SU (2) transformation on q is written 



(2) 



where 



rotation about an axis lying in the spin x-y plane. Defin- 
ing the unit vector Qi — x Qi) = (— sin cos 0;, 0) 
we have 



exp[z7 



cos (hOi) + i sin (i0i)wi • a, 



(4) 



where cos 6i = 2" • fi^ . 

We then find the following Hamiltonian in the new spin 
coordinate system 



lf = ^[-(f^cLM^c,^-Hh.c.) 



(ij) 



+r^iSrS^Q%-\n.n,) 



(5) 



where we have introduced AP^ = (f/^ ' '3'"' = 

RZ^R^ , and {R^)ij = cos6Sij -I- (1 — cos6)uJiUjj + 
sinOeijkU^^ is the 5'0(3)-rotation operator induced by 

Until this point, the discussion has been completely 
general. In order to make further progress, we neglect 
spin fiuctuations and make the restriction that a site is 
occupied by at most one spin and that the electron which 
occupies site i has its spin pointing along the local posi- 
tive z-axis, fli. This contraint reduces Eq. (^ to 



off 



E 



[-(t'^c|cj + h.c.) + K'-^n^n.j 



(6) 



(3) 



The specification of fixes U-^ only up to an overall 
rotation about the new local z-axis. Choosing Gi = 
exp[— zia^rJi • a] makes U!^ — GiU-^ also satisfy the 

defining relation, Eq. (H). To fix this remaining degree 
of freedom, we arbitrarily choose to correspond to a 



with = t^^Mi^^, K'i = \ J^'(^i ■ - 1), and Cj = Cji. 
At this point we have an effective model describing "spin- 
less" fermions on a lattice with hopping amplitudes and 
interaction-strengths being functions of position. In this 
way we have automatically solved the constraint of no 
double occupancies at the expense of treating the spins 
as classical variables and ignoring their quantum spin 
fluctuations. 
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Let us now turn to the properties of t'-' . First we 



note that 



(t^^)*. The complex phase of r*-' can- 



not in general be gauged away by a local transformation 
Cj — > e^'^^Cj if the spin configurations are noncoplanar. 
When we defined the local transformation U-^ , we noted 
that it was only specified up to a rotation about the local 
z-axis, U-^ — > ^^Q- "^^^ effect of such a local rotation 
on the effective hopping cannot be distinguished from 
a local gauge transformation c. — > e'^^'c. which does 
not affect the physics. Hence, the set of physically in- 
equivalent choices of C/jj belong to SU{2)/U{1) = S*^ 
and we conclude that in the absence of an external mag- 
netic field, there are two physical degrees of freedom per 
site (or plaquette) which determine the effective coupling 
constants. 

To understand more precisely which portions of the ef- 
fective hopping r'-' are U{1) gauge invariant, we look at 
a particular plaquette of the lattice, consisting of points 
labelled counterclockwise as ro,ri,r2, and and with 
associated spins pointing in the fli directions. The flux 
through the plaquette is given by 



plaquette 



lmln(r''«''iT''i''=^r''=^'^^T''^''"). 



(7) 



Topological arguments show that ^plaquette — 
^^(ilo, i^i, i^2, f^s), half the solid angle enclosed by the 
shortest path on the sphere connecting the vectors {fti}, 
i.e. the flux corresponding to the plaquette is equal to 27r 
times the skyrmion charge represented by the placuiette. 
Generalizing a result of Wen, Wilczek, and Zeea who 
computed the imaginary part of the formula below one 
may show that for a plaquette consisting of exactly three 
sites 



(8) 



Assuming smooth fields 6{r) and 0(r) one finds in the 
continuum limit that the fictitious flux corresponds to a 
Berry gauge field 

= =^sin2^9^(/), 

which is the vector potential due to a magnetic monopole 
of strength — i. The local U{1) degree of freedom, repre- 
sented by the G's, corresponds to a gauge transformation 
of the topological vector potential A^j,. 



The case of the ferromagnet 

Let us now consider the plaquette (fio, f^i, f^2, f^a) in 
spin space, drawn as in Fig. ||, the parallelogram repre- 
senting a patch of the surface of the sphere. To simplify 
the argument further, we restrict ourselves to the case 
fli ■ Qj = cos 6 for all nearest neighbor pairs (ij) in the 



plaquette, so that the relative angle between the spins on 
each side of the square is 0. A straightforward application 
of spherical geometry yields the relation ^ A = (q;+/3) — tt 
where a and (3 are the interior surface angles of the spher- 
ical parallelogram. 
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Z=2(a+p) 



FIG. 1. The close to ferromagnetic spin-configuration is 
shown in spin-space. All sides in this spherical parallelogram 
correspond to an opening angle 0. E denotes the sum of the 
angles within the path 0^1— >2^3— >0. 



However, since 9 is fixed on all sides, we cannot freely 
choose a and f3. Rather, we solve for <i> = ^A in terms of 
a — p and and then derive expressions for the nearest- 
neighbor coupling constants r and K in terms of the 
same variables. Using standard formulas from spherical 
geometry together with Eq. we find 



I T„„ I = cos(i6') 

$ = 2arcsin (tan^(i6') cos(i(Q; 
K = ij(cos6'- 1)," 



m) (9) 



where r„„ is the effective nearest-neighbor hopping am- 



plitude. We note that $ is given by $ 
for small 6. 
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C0s(i(Q;— /3)) 



The case of antiferromagnetic order 

Let us now consider the analogous calculations for a 
spin configuration which is close to antiferromagnetic. 
Let us denote the local antiferromagnetic spin as fli as 
before. Denote by n[ = —ili the antipodal points. Not- 
ing that in the case of a spherical parallelogram, great 
circles which connect the four sides of {Qi} intersect in 
the points {f^^}, we find that the path fio — > il[ fl2 ^ 
f2g — > f2o connecting the antiferromagnetic spins has the 
geometry sketched in Fig. ^. Fig. ^ shows another il- 
lustration of the path taken in spin space when going 
around a plaquette in an antiferromagnetic background. 
By using Eq. ( ^) the angles a have become exterior 
rather than interior angles in the path. This gives us the 
relations 



I '^nn 

K 



sin(i6') 

n-{a-(3) 

-jj{cose + 



1) 



(10) 
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where in this case cos 9 = — ilj • ^i^i- Comparing to 
Eq. ( ^ we see that the connection between 9 and $ has 
disappeared. 




3,-' 
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I=47t+2(|3-a) 

FIG. 2. The close to antiferromagnetic spin-configuration 
is shown in spin-space. The left part corresponds to the points 
Qi, while the right part corresponds to the antipodal points 
It is seen how the Q-angles become exterior rather than 
interior leading to a completely different E as compared to 
the ferromagnetic case. 



By extending the same argument to the adjacent cell, 
we reproduce the arguments, however all paths are tra- 
versed in the opposite sense and the flux becomes neg- 
ative compared to the first cell . The flux is therefore 
staggered, with a ± sign depending on the sign of the 
sublattice associated with the plaquette. 




FIG. 3. The path taken in spin-space when going around 
a plaquette ^ I' — > 2 3' — > in an antiferromagnetic 
background. All paths follow great circles on the sphere. 



In contrast to the case of the ferromagnet, we find that 
we can change the effective staggered flux for the electron 
gas, Eq. (|^), from — tt to tt without affecting the magni- 
tude of the local values of the effective nearest neighbor 
coupling constants. In fact, for the most symmetric case 
a = /?, we find that an antiferromagnetic skyrmion gener- 
ates an effective staggered flux of exactly tt per plaquette. 
Due to the curvature of the sphere, such a texture cannot 
in general be extended to the plane. 



Next-nearest neighbor hopping 

When we include second neighbor hopping, the situa- 
tion becomes much more complicated for the antiferro- 
magnet. Each square plaquette has four gauge invariant 
fluxes, corresponding to each triangle defined by the re- 
moval of one vertex from the four corners of the square. 
When the spins are aligned close to ferromagnetically, 
the shortest path in spin space connecting the two di- 
agonal spins lies wholly within the region defined by the 
corner spins. Hence the flux through each triangular sub- 
plaquette is very closely proportional to the area of each 
subtriangle in real space, and the sum of flux through 
a pair of subtriangles that cover the square is precisely 
equal to the flux through the entire plaquette. As a con- 
sequence, the effective phase of each of the interactions is 
extremely well approximated by spreading a constant ef- 
fective magnetic fleld corresponding to the local skyrmion 
density throughout the entire real space plaquette and as- 
signing the effective phase on the links by a conventional 
choice of gaugeO. 

In the case of a conflguration close to antiferromag- 
netic, the four subtriangles cover the entire sphere in spin 
space [see Fig. ^ . Taking into account the orientation of 
each of the bounding paths, we find the following relation 
between the flux through each of the triangles for both 
the ferromagnet and antiferromagnet: 

- $(17of^i'f^3') - ^{hi'h2^3') = 27rn, (11) 

where n ^ for the ferromagnet and n = ±1 for the an- 
tiferromagnet. Note that Eq. ( ^l|) is valid in general, i.e. 
it does not rely on the assumption of a spherical paral- 
lelogram. It is also easy to see that it is valid in the pres- 
ence of an external electromagnetic flux. If we assume 
a spherical parallelogram we can also note from Fig. ^ 
that $(f2of2i/f]3/) = $(f2i/02l^3') = 5$. It therefore fol- 
lows that ^{Qi>Q2^o) = ^{flo^2^3') = |$ + n7r. Hence, 
for the antiferromagnet, links connecting sites within the 
same sublattice pick up a phase tt if they belong to sub- 
lattice,A, and zero if they belong to sublattice B, or vice 
versa.l3 

In Fig. ^ is drawn a 3 x 3 subset of a 2D antiferro- 
magnetic lattice. The value of the local spin generated 
phase Imln(M-[^) = is indicated on each link. Ar- 
rows indicate the direction i j. The gauge is chosen 



so that 



0. The phase 



_ (-1) 



$. Along 



the diagonals, the phase of the coupling constant is: 

- m + i-iy] 

^ r,r-\-y-\-x 2.^ ^ ' ^ 

(\> = f [1 - (-1)1. 



(12) 



The local gauge choice does not take account of the 
curvature of the spin space and in general, the sphere 
cannot be covered with parallelograms so the local choice 
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of gauge cannot be extended to the whole sphere. Never- 
theless, in the limit — > the above formula is expected 
to be valid. 




FIG. 4. A particular choice of gauge for the staggered 
flux-phase (±4> ) with next-nearest neighbor hopping is in- 
dicated through the phases associated with each link in the 
lattice. 

An intriguing special choice for 9 where it is possible 
to cover the sphere by parallelograms is cos 9 = —1/3. 
In this case, the spins reside on the vertices of a regular 
tetrahedron, and it is also possible to extend this map- 
ping exactly over a square lattice. Considering that going 
around a plaquette in the real space lattice, corresponds 
to encircling two of the faces of the tetrahedron, it is ob- 
vious that this loop will cover half of the spin-sphere and 
hence the "tetrahedral phase" generates a staggered flux 
±7r. 

Putting all this together, we find that given origi- 
nal nearest and next nearest neighbor hopping 
(tnnn), the cffectivc coupling constants for the remain- 
ing fermionic degrees of freedom are determined by the 
two free parameters 9 and which control the unstag- 
gered and staggered fiuxes, respectively. In the general 
case with an external electromagnetic flux we have 2 + 1 
physical degrees of freedom per plaquette. This results in 
three independent fluxes per square, the only restriction 
being given by Eq. (|lll). Using standard formulas from 
spherical geometry, we find that with a suitable choice 
of gauge, the effective nearest and next nearest neighbor 
couplings are given by: 



CONCLUSION 

The main conclusion from this paper is that in con- 
trast to the ferromagnet, antiferroiHagnetic spin textures 
naturally generate a staggered fluxQ on the basic nearest 
neighbor plaquettes. However, in contrast to case of the 
ferromagnet, the amount of flux per plaquette is weakly 
related to the magnitude of the effective nearest neigh- 
bor coupling. Incorporating second neighbor interactions 
reveals the existence of an additional independent effec- 
tive U{1) fictitious flux on plaquettes of second neigh- 
bor sites. Since the flrst neighbor coupling is severely 
suppressed in a Neel antiferromagnet, it is quite possi- 
ble that the effective second neighbor interactions could 
be important; in the limit of weakly coupled effective 
nearest neighbors and strongly coupled second neighbors 
textures such as appear iiL.the two layer quantum Hall 
problem could be favored. El However, even in this case 
the magnetic field energy generated byrf taggered orbital 
currents could supprjess these textures.tl 

Using Eq. (|l3|) well3 have numerically calculated ener- 
gies of various uniform spin textures in the Hartree-Fock 
approximation and for nonuniform stripe textures in the 
Hartree approximation, with the goal of understanding if 
noncoplanar spin textures which generate fictitious stag- 
gered flux have lower energy than the coplanar "^.iral". 
phases studied previously by a number of authors.l30'E£l 
Preliminary results are consistent with those that have 
been obtained previously by more elaborate means; spi- 
ral phases appear to have the lowest energy for the uni- 
form phases but are thermodynamically unstable against 
phase separation. The noncoplanar phases all have stag- 
gered fictitious rather than the unstaggered flux which is 
known to lower the kinetic energy in an electron gas and 
there does not appear to be any clear association between 
topological charge and doping, at least in a tight bipding 
model dominated by nearest neighbor interactions.EJ 

By neglecting quantum spin fluctuations as we have 
done, our numerical simulations of the nonuniform spin 
textures suffer from the same problem as others who sim- 
ilarly ignore spin fluctuations. The ordinary Neel anti- 
ferromagnet becomes insulating, and this, together with 
the tendency towards phase separation apparently does 
not favor noncoplanar spin textures with any nontrivial 
topological properties. It remains to be seen if properly 
accounting for quantum spin fluctuations can change this 
result. 



'^nn ^nnSin(2^) 



„i(-l)'-$/2 

(13) 
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